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' Abstract. The status of stopping in heavy ion reactions is reviewed by comparing 

, available data in pp, pA and AA systems over the energy regime ~ 2.5 — 130 

' GeV. The data consist of average rapidity losses, anti-proton over proton ratios, and 

\^ , net-baryon values at mid-rapidity. The overall features of nuclear stopping are reason- 

• bly well described by simple extrapolations of pp, and pA collisions to AA. 
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1. Introduction 

The major goal for relativistic heavy ion reactions is to form hot and dense nuclear matter, 
. , to study its properties, and to come to a better understanding of non-pertubative QCD. 

' The heavy ion experiments at RHIC have this as their main purpose. It is important to 

H I understand nuclear stopping in detail since it is a requisite for the formation of such systems 

in describing the conversion of the initial kinetic energy into matter excitation at mid- 
rapidity. A description in detail of how this transport takes place is a necessary ingredient 
in our overall understanding of the reactions, and in the expectations to form and study the 
properties of QGP. 

The importance of stopping was recognized early [ |lj ^ with a good conceptual un- 
derstanding of pA and AA stopping. This also laid the foundation for estimations of what 
energy densities might be reached in heavy ion induced reactions. Experimental data have 
come at a slow but steady pace as the accelerators and experiments have been build and 
data analyzed. Some years back I surveyed the available data with Ole Hansen [ ^, and 
now is a good time to revisit the subject, since more detailed data has emerged at lower en- 
ergies, not only in AA collissions but also in pp and pA. In addition, the first run at RHIC 
in 2000 has given us important clues about the de facto situation on stopping at extreme 
ultra relativistic energies. In this paper I will review some recent experimental data on stop- 
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ping coming from pp, pA at lower energies, as well as data extending our knowledge about 
AA. I will put the data into perspective by comparing rapidity losses with estimates from 
the multichain model [ ^ for the energy dependence. Though the first data from RHIC 
did not measure baryon rapidity distributions and thus enable us to quantify the amount 
of stopping, the particle production at mid-rapidity and the anti-proton to proton ratios are 
indicator of stopping and is discussed in the last section. 



2. pp and pA reactions 

It has often been conjectured that a thorough understanding of pp and pA reactions is 
needed to gauge the significance of results from AA, whether results can be interpreted as 
merely simple superpositions of more elementary collisions or if other many-body effects 
play a role. Such superpositions or extrapolations clearly rely on high quality data from 
pp and pA. The data that formed the basis for the early understanding were in many ways 
limited, in large part due to lack of statistics and lack of centrality event selection. I will not 
review these here, but merely point to the references given in the earlier review [ |l], |[ ^ ^]. 
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Fig. 1. dN/dx/ distribution for pp, and pPb from preliminary NA49 data [ 



In recent years significant data have emerged from experiments at AGS and CERN. In 
particular the data from E910, E941 at the AGS, and NA49 at the SPS ai-e most significant 
in terms of rapidity coverage and centrality selection. The data have yet to be published, 
but preliminary results have been presented at the Quark Matter 2001 conference in Stony 
Brook [ § 0, and at this conference [ In general the data have been selected on 
centrality by counting 'grey' tracks i.e. slow protons emitted from the target. The 'grey' 
tracks have been shown by model calculations to be related to the number of collisions that 
the incoming proton suffers in the nucleus. The data show that after the first collision the 
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system has forgotten the isospin content of the entrance channel as manifested by the ob- 
servations that spectra, and rapidity distributions of pions, protons and neutrons are similar 
[ ^ when gated on number of collisions greater than one. 

The NA49 experiment has made further progress by isolating the behavior of the 
stopped projectile by nucleus. This is done by subtracting the component from the target 
measuring the protons in 7I-A reaction. The pion has no baryon content, and the produced 
protons correspond to the target contribution in pA reactions. The preliminary dN/dxp 
distribution are shown in Fig. The Feynman xp is defined in the interval -1 to 1 as 
Pi(c.m.)/Pi(c.m.)„,„;f The dN/dxp for pp is almost constant, as has been known for a long 
time. The pPb data are shown for two values of the mean number v of collisions, v = 3.1 
corrsponding to mean bias pA and v = 6.3 to central collisions. The central collision events 
exhibit complete stopping with the peak of the distribution at Xp =0. These detailed mea- 
surements confirm the previous expectations of Busza and Goldhaber [ |l[]. NA49 also 
measured distributions of neutrons emitted in pp and pA and found the distributions after 
one collision to be nearly identical to those of stopped protons. 




Fig. 2. Pion and proton dN/dy distributions in NiNi reactions at 1 .93AGeV from the FOPI 
Collaboration [|l^. 



3. AA Collisions 

An example of the data obtained since the previous review is shown in Fig. ^ the rapidity 
distributions of pions and protons in NiNi reactions at 1.93 A GeV [ |lO||. The data are for 
three very central collision bins. The most central 100 mb cross section bin show a slight 
enhanchement both of pion production and of protons, indicative of a slight increase in 
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stopping. The data are compared to the expectation of a pure isotropic Boltzman distri- 
bution with a temparature of T=125 MeV. Already at this low energy we do not obeserve 
complete stopping in the sense of isotropic decay of a fireball, but see effects of incomplete 
stopping or of flow in terms of longitudinal expansion. 

The E917 experiment has extended the earlier measurement of E866 [ [ll|] at 11.7 
AGeV/c with measurements at 6, 8 AGeV [ |l2|]. The proton rapidity distributions are dis- 
played in the left panel of Fig. |[ Again the data are compared to expectation of isotropic 
decay, and incomplete stopping is observed. The data were analyzed in a novel way, at- 
tributing the overall distributions as being a sum of two Gaussians representing an incom- 
plete stopping of projectile and target protons. Such decomposition attempts to determine 
the projectile stopping as expressed by the average rapidity loss. The values extracted by 
this method are similar to those obtained and described in the next paragraph. 




-2-1012-2-1012 



Fig. 3. The left panel shows proton rapidity distribution at AGS energies from E917 [ [T^ ], 
and the right panels the extracted temperatures from Boltzman fits to the data. 



For symmetric systems the average rapidity loss of protons/baryons are defined as 

JdydN/dy{yp-y) 
^''^ = JdydN/dy 

The integral is taken from mid-rapidity (y= 0) to the beam rapidity yp and is an approxima- 
tion to the actual projectile projectile, since it is not possible to distinguish between target 
and projectile protons when protons pile up at mid-rapidity. Figure ^ symmetrizes the re- 
sults of AA average rapidity loss byp relative to the beam rapidity versus the kinetic energy 
per nucleon. With the additional data available it is apparent that both the heavier systems 
and the lighter systems both have constant relative rapidity losses with the heavier slightly 
higher with an average value of 0.32. Such higher value is of course not unexpected due to 
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Fig. 4. Average rapidity relative to the beam rapidity vs. laboratory energy. The experi- 
mental data are for light systems SS SiAl, and NiNi, and heavy AuAu and PbPb, from the 
FOPI, E802, E866, E917 and NA49 experiment. The NA49 data ai-e from Ref.[ [l3]l. The 
reference for the other data are given in the text. 



the increase in thickness of the nuclei involved. Such plot is also deceiving in displaying 
no energy dependence. The absolute value does increase with energy, and so the amount 
of available energy for particle production. The expectation of Sy^ from a isotropic Boltz- 
man distribution with a temperature of 125 MeV is shown on the figure rising from about 
0.30 at the lowest energy, again illustrating that such isotropic stopping is not observed at 
relativistic energies. 

4. Multichain Model 

The multichain model was introduced in Ref.[ |2| to study stopping in pA and AA reac- 
tions. This model is used to illustrate how well the overall features of nuclear stopping is 
described by such simple approaches. Similar analysis could also be carried out using one 
of the many event generator available. These models have quite similar assumptions. The 
basic assumption and points are 

• The nuclear geometry is seperated from the dynamics. 

• Transverse distribution are independent of E,A, and p/,. 

• Glauber or wounded nucleon model is used to calculate collision probabilities. 

• Multiple collision dynamic scales with Feynman x. 

• Independent projectile and target fragmentation. 
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The calculations presented here utilize equations presented in section 2 and 5 of Ref. [ 
with a constant fragmentation function for pp, and and a value of a=1.3 at the lower 
energies and 2.7 at the higher. Such value was derived in [ ^ for AGS and SPS energies, 
respectively. Since the calculation here is done down to low values of ^/Inn , both target 
and projectile contribution were included, and at cutoff was introduced in the allowed final 
XF to require open inelastic channels. 




Fig. 5. Calculated rapidity density distributions fo'if 4 cm. energies using the multichain 
model. 



This relatively simple model allows us to make predictions of the rapidity loss vs cm. 
energy. The results from this calculation, assuming a standard Wood-Saxon distribtuion 
of the nuclear densities, and an inelastic nucleon-nucleon cross section of 30mb at the 
lower energies and 40 mb at RHIC energies are displayed in Fig. |^. The 4 panels show 
un-normalized rapidity distributions for the projectile and target baryons as the dashed 
lines, while the total distributions are given as solid lines at four cm. energies from SIS to 
RHIC. The evolution from overlap between projectile and target, and nearly full stopping 
to partial transperency is clearly observed. The average rapidity loss byp is calculated 
from the distributions and the plotted relative to the beam rapidity in Fig. ^ versus the 
cm. energy. This schematic model does in fact quite accurately describe the observed near 
constant value in the lower energy data shown in Fig.^. Similar features are also predicted 
from e.g. cascade code, but the present model analysis demonstrates that the main features 
are described starting with the elementary pp energy loss mechanism, and folding with a 
geometric description of heavy ion reactions using a Glauber or Wounded Nucleon Model. 

5. Early Results from RHIC 

Most of the early results from RHIC have been focused at mid-rapidity, determining charged 
particle densities, and for identified particles looking at ratios to minimize the systematic 
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Fig. 6. Predicted energy dependence for the average rapidity loss in the multichain model. 
The dashed line shows the near constant value observed in AA collisions at lower energies. 



errors. Thus measurement of net-baryon number densities at mid-rapidity (though prelimi- 
nary numbers have been quoted by STAR), nor a determination of the rapidity distributions 
of net-baryons from which rapidity losses can be determined has been made. Even so 
these measurements shed light on stopping at ^/snn =130 GeV by investigating the ratio 
N(p)/N(p) at mid-rapidity and its centrality and pt dependence. So far 3 experiments have 
reported final values of N(p)/N(p) (see Refs [ [T^ [l^). The values all agree and are 
in the range of 0.57 — 0.65 for central collisions neary = 0. It is of interest to compare this 
value with heavy ion data from lower energies. Data exist for Au-Au collisions at the AGS 
[ |l7[l, near threshold for p production, and at SPS [ |l8|]. Data at higher energies have been 
measured in pp reactions at the ISR [ |l9[ These later data cannot directly compared 
with the heavy ion data, but require an iso-spin correction. A simple estimate for such 
correction is described in the appendix. The available data are displayed in Fig.^. The ISR 
data are displayed both with the measured values (open symbols) as well as with the isospin 
corrected (filled symbols). An overall smooth dependence with ^/Inn is observed with 
the RHIC heavy ion data being the largest. The values observed at ^/snn =130GeV are 
also close to an extrapolation of the N(p)/N(p) ratios from ISR [ pT| ]. 

BRAHMS has measured N(p)/N(p) not only at y=0, but also at two larger rapidity 
values [ |l5|], which shed further light on the issue of stopping. The rapidity dependence of 
the N(p)/N(p) ratio is shown in Fig. ^ The present systematics show that the ratios fall off 
more rapidly over two units of rapidity than the nucleus-nucleus results at lower y^s^N, but 
are very similar to the pH-p result at roughly half the CM energy [ 20 1 . 

Also in Fig. ^is compared the measured ratios to calculations using the HIJING model 
[ Hi, the FRITIOF 7.02 [ ||] string model and the UrQMD cascade model [ ^ using 
the same cuts on centrality as those applied in the data analysis. The FRITIOF model 
reproduces the N(p)/N(p) ratios quite well, while it overpredicts by « 30% the charged 
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Fig. 7. Energy dependence ofp/p ratio from AA and pp data. 



particle yield at r| w [ |2^]. This is mainly related to a fairly large degree of stopping 
predicted by the model. On the other hand the HIJING model descibes the overall charge 
particle yields at r| w while it fails in describing rapidity dependence of the anti-proton 
to proton ratio. This feature of the model is in turn related to a small stopping of the 
projectile baryons. The UrQMD model, which is not a partonic model, underpredicts the 
ratios by almost a factor of two. These models exemplify the present understanding of 
heavy ion collisions at this new energy regime; none of the models offers a consistent 
description of the observed features, and the rapidity dependence serves as a benchmark 
for the description of stopping and anti-baryon production. 

6. Conclusions 

The issue of stopping in heavy reactions has been studied in details in the past fifteen years 
and systematic heavy ion data on this exists from SIS to SPS energies, and will emerge 
from the data at RHIC within the next year The rapidity losses have been measured as 
functions of projectile, centrality, and energy. The existence of new high quality data from 
pp and pA will enable us to see if the AA stopping is merely simple extensions of pp and 
pA, or if new features due to the many-body overlaps manifest themselves in the details. 
The data reviewed here call for a careful theoretical analysis, and not a phenomenological 
approach. 
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Fig. 8. Comparison of the measured N(p)/N(p) ratios to model predictions (BRAHMS 
data). The data shown are for — 40% central events. The three model calculations (HI- 
JING, FRITIOF, and UrQMD) are shown for comparison. See text for details. 



7. Appendix 

An isospin correction is needed when comparing N(p)/N(p) ratios from pp with data from 
AA. The correction arises because the p multiplicity from pp-^ pX is 0.6, while for an 
isospin symmetric Nucleon-Nucleon (NN) system we expect NN^ pX to have a multi- 
plicity of 1.0. The first order correction is estimated as follows. The proton cross section 
as coming from two components a) the direct stopping Os and b) a pair-production compo- 
nent Op giving rise to an equal cross section of protons and p's. The pp reaction measure 
Op I {<3p+(3s) , while for the comparison with AA we want to estimate the same quantity 
with isospin weighted yields. This is achieved by replacing the Os with a increased yield of 
+<ys/Q.6 In that case the isospin weighted ratios On / (Op+Gs) is given by Gp I (a/>H-as/0.6) 
with yields from the pp data. The expected ratio for NN collisions can then be written as 
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A^(N)/A^(N) = A^(p)/A^(p)/(5/3 - 2/3 * A^(p/A^(p)) 



(2) 



This has been used to correct the N(p)/N(p) ratios from ISR in Fig.^. 
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